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Abstract 


This note records the solution or partial progress made on several open 
problems related to poly-Bernouilli (pB) numbers, posed by Knuth in [1]. 


1 Problem 1 


Is there a combinatorial proof for the following formula? 
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The left part of the formula is the number of parades with m+ 1 girls and n 
boys starting with a girl. Those are all the parades of the form agi. We put 
the youngest girl g, in our parade. For a we choose m — k of the m remaining 
girls and n — 1 of the n boys. The remaining & girls, | boys come after the 
youngest girl and make up @. Now a has to be a parade starting with a girl or 
be empty, because our full parade has to start with a girl. It has to end with a 
boy or be empty, because otherwise g; is preceded by an older girl, which is not 
valid in a parade. So the number of possible a is Dm—zn—1 by equation (17.1) 
in [1]. For 6 any valid parade will do, so there are B;,; possiblilities. And there 
are ie) (7) ways to divide the girls and boys between a@ and (. 

Using variations of the above combinatorial argument, we can dedude new 
identities between Byin, Cm n and Dm». For example, by using gm, the oldest 
girl instead of gi, we find that 


C2 = BPD (‘") G) (Ce_ait[b = 1 = 0])(Cm—an—1-1+[m—k = n—I-1 = 0])) 


If we start with 6), the youngest boy instead of g1, we find (using [17.2] from 


[1]): 
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And if we start with 6, we find 


m+lin-1 
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Many more identities can be derives using similar constructions on other 
classes from section 17 in [1]. 


2 Problem 2 


All mxn matrices for m < 5,n < 5 were enumerated using a computer program. 
Then for all (?) matrix pairs {¢, } the program counted how many of the mxn 
matrices were {¢, w}-free. There were 25 pairs that resulted in the pB-numbers. 
(See table 1). 

Bijections for most of these {¢, W}-pairs were already discussed in [1]. One 
slighty new variation are the pairs in the min/max-closed section where one of 
the matrices is € a): The pair { € i) ({ i) \ corresponds with the 
*max-min’-relations, which satisfy the condition 


xi, yj and xy ~ yj implies Umax{i,i/} —~ Ymin{j,j’}- 


And similarly, { (; a (i i) } corresponds with the ’min-max’-relations, 


which satisfy 
Ui Yi and xj ~ Yi implies Lmin{i,i/} —~ Ymax{j,j’}- 
We can construct bijections to the [-free matrices by repeatedly changing 


2x2 submatrices of the form & : to ({ ) (respectively changing & 7 


3 Problem 5 


Investigate the behaviour of the nonuniformity ratio maxo<r<nPmn(r)/ mino<r<nPmn(T) 
when m and n are large. 
Experiments seem to indicate that, for a given n, if m is large enough (~ m > 
n/2), the maximum Pyn(r) is reached when r = 3, the minimum is reached 
when r = 0. Can this behaviour be explained? 


4 Problem 9 


This problem asks to investigate the number of lonesum / x m x n tensors Lj mn 
and also the weakly lonesome | x m x n tensors Limjn- 


In bijection with 


Min/max-closed relations 
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Table 1: The 25 pairs {¢,} that result in pB-numbers. The right pair on each 


line is equal to the left one with the ones and zeros swapped. 


For the 2 x 2 x n weakly lonesum tensors, we only have to make sure that 
A= ; : and B= : i) don’t occur together in the same tensor. Let’s 
call s,(U) the number of tensors of size 2 x 2 x n that don’t contain A and B, 
$,(A) the ones that contain A and s,(B) the ones that contain B. This gives 


the following recurrences: 
so(U) = 1, s0(A) = s0(B) = 0 
sn41(U) = 148,(U) 
8n41(A) = 8n(U) + 158,,(A) 
$n41(B) = 8,(U) + 15s,(B) 
fn = 1+ 15s,(A) + 15s,(B) + 16s,,(U) 


Solving these recurrences gives the generating function: 


13z-1 


Goon(z) = (1 — 15z)(1 — 14z) 


This results in the following formula: 
Loon = 2x15" -—14" 


Using a similar method, we can also find the generating function and explicit 
formula for the 2 x 2 x n lonesum tensors: 


1 — 212? — 2823 


Loon = 8x 5"—4x 4" —4 x 3% 42” 


5 Problem 10 


This problem asks to investigate the number of max-closed | x m x n tensors. 
We will compute the generating function for the 2 x 2 x n case. 

If we start with a 2 x 2x 1 tensor and try to extend it to size 2 x 2 x 2, there 
are 14 possible 2 x 2 x 1 tensors that can occur to the right as a successor. Let’s 


write abcd for the tensor : i) Let’s represent the successors S(t) by a bit 


string of length 14, with the bit positions corresponding to the tensors as follows: 


1111, 1101, 1100, 1011, 1010, 1001, 1000, 0111, 0101, 0100, 0011, 0010, 0001, 0000 


There are only 7 possibilities, shown in table 2. 


Tensor t | Possible successors S(t) | Count c(t) 
*000 11 1111 1111 1111 14 
*001 11 0101 0110 1011 9 
*010 10 0110 0110 1111 9 
*O011 10 0100 0110 1011 7 
*100 11 1000 0111 1011 9 
*101 11 0000 0110 1011 7 
*111 10 0000 0110 1011 6 


Table 2: List of possible successors for each 2 x 2 x 1 tensor. 


If we have a tensor t,t ...tn_1 made of 2 x 2 x 1 tensors t;, we can extend it 
with a tensor t, where t,€ {) S(t). 
i=1..n-1 
Let call s,(t) the number of max-closed 2 x 2 x n tensors that have the set 
S(t) as a successor. Then we can setup recurrences for the s,,(t), and based on 


those count the number of 2 x 2 x n tensors fr. 
80(*000) = 1, so(t) = 0,¢ # «000 


Sn.1(*000) = 25, (#000) 
(001) = 25, (*000) + 35n (+001) 
1(*010 ) + 38, (*010) 
$n41(*011) = 28,,(*000) + 28,,(*001) + 38,,(*010) + 4s,,(*011) 
( ) (100) 
(100) 4 


Snt 


Sn4 


$n41(*100) = 28,(*000) + 3s,(*100 
8n41(*101) = 28,,(*000) + 25, (*001) + 38,,(*100) + 4s, (*101) 
Sn41(*111) = 28p (+000)+28,, (*001)+3sy (+010)+38,, (*011)4+-35, (+100)+3s,, (*101)+65, (*111) 


fr =14+ > Sn-1(t)e(t) 


Solving these recurrences gives the generating function 


Foon(z) = 


(1 — 6z)(1 — 4z)(1 — 3z)(1 — 2z) 


From this generating function, we can find a general formula for fo2n: 


fot = 35K 6 = 24 x a ek Bt 


Using a similar method, it is possible to find Fy, and fimn, for a few other 
Ixmxn: 


(1 — z)(1 + 2z — 3092? + 2846z3 — 920824 + 8516z° + 34082°) 
(1 — 10z)(1 — 8z)(1 — 7z)(1 — 6z)(1 — 5z)(1 — 4z)(1 — 3z)(1 — 22) 


Fo3n(Z) _ 


foan(z) = 1547156x10"—* 


203 105 74 
50 x8"—44x7"4 . x67+432x57 147x4” 5x3" x2” 


1 + 42z — 42112? + 10891423 — 129118724 + 6343786z° 
+138616072° — 3444887342" + 18324013902 
(=z) — 455000512829 + 5074183440z!° — 1696749984211 
(1 — 152)(1 — 132)(1 — 12z)(1 — 11z)(1 — 10z)\(1 — 92) 
(1 — 8z)(1 — 7z)(1 — 6z)(1 — 5z)(1 — 4z)(1 — 3z)(1 — 2z) 


Poan(z) == 


From the above, we can conjecture that Fo,,,, has the form 


(1 — z)(some polynomial in z of order() — 4) 


Tkat.,...."2)-2,("2) 1— kz 


1+ 90z — 101572? + 3476702? — 557799724 + 340443142° 
+2725072812z° — 7049450706z" + 6269000036028 — 307963368792z° 


+ 869031680800z19 — 12858128750722!! + 701635229056z!? — 1116887393282!3 
F33n(z) = (1 z) 


(1 — 202)(1 — 17z)(1 — 14z)(1 — 132)(1 — 122)(1 — 11z)(1 — 102) 
(1 — 9z)(1 — 8z)(1 — 7z)(1 — 6z)(1 — 5z)(1 — 4z)(1 — 3z)(1 — 2z) 


6 Future work 


Here are some additional open problems (some of which may be answered in a 
future version of this note): 


1. Are there any other interesting numbers and corresponding bijections hid- 
ing in the 120 — 25 = 95 pairs of matrices that don’t give pB-numbers in 
Problem 2? 


We can compute the number X,,,, of {¢, w}-free matrices for all 95 pairs {¢, W}, 
and then group them in classes with the same X77, for all m and n. Then there 
are (1, 5, 2, 1) classes of size (1, 4, 8, 16) where Xm» = Xn,m, and (2, 6, 2, 0) 
classes where Xmn # Xn,m- 

The symmetric class of size 16 seems a good starting point for further inves- 
tigations. 


2. Is there an algorithm to generate all parades of m girls and n boys by 
swapping 2 persons at each step? 
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